Josephson effect through a quantum dot array 



O 

o 



(N 



u 

Oh' 



o 



X 



F.S. Bergeret, A. Levy Yeyati, A. Martm-Rodero 
Departamento de Fisica Teorica de la Materia Condensada C-V, 
Umversidad Autonoma de Madrid, E-28049 Madrid, Spain 

We analyze the ground state properties of an array of quantum dots connected in series between 
superconducting electrodes. This system is represented by a finite Hubbard chain coupled at both 
ends to BCS superconductors. The ground state is obtained using the Lanczos algorithm within a low 
energy theory in which the bulk superconductors are replaced by effective local pairing potentials. 
We study the conditions for the inversion of the sign of the Josephson coupling (vr -junction behavior) 
as a function of the model parameters. Results are presented in the form of phase diagrams which 
provide a direct overall view of the general trends as the size of the system is increased, exhibiting 
a strong even-odd effect. The analysis of the spin-spin correlation functions and local charges give 
further insight into the nature of the ground state and how it is transformed by the presence of 
superconductivity in the leads. Finally we study the scaling of the Josephson current with the 
system size and relate these results with previous calculations of Josephson transport through a 
Luttinger liquid. 
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Advances in the fabrication of hybrid nanostructures consisting of quantum dots, molecules or carbon nanotubes 
attached to conducting electrodes have lead to an increasing interest in the electronic and transport properties through 
'' such low dimensional devices. In most of these systems electronic correlations play an important role leading for 
C/3 , instance to the Kondo effect [H, H, Q. In particular, when the leads are superconductors a natural question that arises 
is whether a supercurrent may flow through the structure. The interplay between Josephson and correlation effects 
have received considerable attention in the last years ^ft, ^5„i|, 7, 8, 9, 10, 1 Ij, [13, • ^ paradigmatic example of this 
interplay is the prediction of a transition to a 7r-junction behavior as a function of some relevant system parameters 
[l4| . This prediction has been recently confirmed experimentally both in carbon nanotubes and semiconducting 
nanowires [3, [1^1 . In spite of all the recent theoretical activity in this field there are still many issues to be explored. 
Q While most theoretical works have concentrated in analyzing the Josephson effect through a single level Anderson 
O . model coupled to superconducting electrodes, the analysis of more complex models taking into account the spatial 
extension of the correlated region is still lacking. In this respect some works addressing the double quantum dot 
situation have appeared recently [ll|, . A natural extension of the Anderson model to describe a correlated region 
^ ' with a finite size is provided by the Hubbard model. Such a model could describe different physical situations like 
' quantum dot arrays, molecular junctions or nanowires. In the present work we will take the linear Hubbard chain as 
0^ ' a prototypical example. Some aspects of this model have been recently discussed in Ref. fl^. On the other hand, 
studies of Luttinger liquids in contact with superconductors [19, 20] provide predictions on the Josephson effect which 
, the model should obey in certain limits. As we show in this work the superconductor-chain-superconductor system 
t — ' reveals several novel correlation effects which are not found in the strictly zero dimensional Anderson model and which 
' are associated with the finite spatial extension of the system. 

The paper is organized as follows: In section II we introduce the model and discuss the approximation scheme used 
to obtain the ground state properties of the system. We will argue that these properties can be adequately described 
by a low energy theory in which the bulk superconducting electrodes are replaced by local pairing potentials. In 
order to test the validity of this approximation we consider in section III the single site case (single level Anderson 
model) and compare the results with the ones obtained using other methods like Hartree-Fock and finite U slave boson 
' mean field. In section IV we present the phase diagrams illustrating the regions in parameter space where 7r-junction 
behavior is found for different number of sites in the array. Section V is devoted to analyze the spin-spin correlations 
and the spatial distribution of the spin excess which characterize the 7r-junction state. Such analysis is useful for the 
full characterization of the interplay between Kondo, antiferromagnetic and pairing correlations. Finally in section VI 
we present results for the current-phase relation and its scaling behavior as the number of sites is increased. By means 
of these results we discuss the connection with previous studies for Luttinger liquids connected to superconducting 
leads The paper is closed with some concluding remarks. 
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II. MODEL AND METHOD 

In the present work we shall assume that the system can be described by the following Hamiltonian 

H = Hl + Hr + Hchain + Ht (1) 

where H]^(^fj^ corresponds to the left (right) superconducting lead, characterized by order parameters A^(j^')ea;p(i0/^(^)) 
and Hchain denotes the finite Hubbard chain Hamiltonian with N sites 

N N N 

Hchain = i^(c|^Cj+lo- + h.c.) + U'^fli^flii + e'^fii^ (2) 
i—l i—1 i—1 

where the field operators creates a particle in site i with spin cr, t is the hopping parameter within the chain, U 
is the on site electron interaction and e is the site energy. The coupling between the chain and the leads is described 
by the term Ht- The determination of the ground state properties of this model is a formidable task which requires 
some approximation scheme. Several theoretical approaches have been used in the single site case {N = 1), including 
numerical renormalization group 0], quantum Monte Carlo [8], finite order perturbation theory Hartree-Fock Q 
or slave-boson mean-field [l3| ■ 

A great simplification can be made by integrating out the electronic degrees of freedom of the superconducting 
leads in the way it was proposed in Ref. [l^. This procedure leads to an effective low energy theory in which each 
superconductor is replaced by a single site with an effective pairing potential A/^(fl;). Also the hopping terms between 
the leads and the chain are replaced by effective parameters tj^^Ry Thus, the terms Hj^^r^ and Ht take the form 

Hl(r) = Ai(^)e^'^-<«)c[-f4^ + /i.c. 

Ht = {iLc\„ci„ + H.c.) + {Irc^j^^crc + h.c.) (3) 

We shall refer to this simplified description of the leads as the zero bandwidth model (ZBWM). In principle the 
effective parameters in this approach have to be determined from the bare parameters by means of a self-consistency 
condition. As discussed by Affleck et al. [l9| using a renormalization group analysis the effective parameters, with 
positive or negative interactions, flow to the condition of zero or perfect Andreev reflection respectively in the limit 
of an infinite chain. However, for the finite chain case these arguments are no longer valid. We shall adopt here the 
simplifying assumption that A^^^ — t^^R — t without an attempt to fine tune them within the range of parameters 
considered. This is a reasonable choice as far as one is interested in the qualitative trends rather than in the 
detailed quantitative results. As we will show in the next section, taking the one site case as a test case, the ZBWM 
yields a phase diagram in good agreement with known results from other approaches that take into account the full 
spectrum of the superconducting leads. One should mention that the Hamiltonian in the ZBWM couples in principle 
configurations corresponding to the grand canonical ensemble. This number can be substantially reduced by 
projecting into subspaces with well defined spin. In order to reach the maximum possible system size we use the 
Lanczos algorithm for the determination of the system ground state . This allows us to obtain results for systems 
up to 10 sites in the chain. 



III. THE SINGLE LEVEL ANDERSON MODEL AS A TEST CASE 



In this section we shall consider the single level case in order to test the validity of the ZBWM for the superconducting 
leads. Previous studies have shown that this system undergoes a quantum phase transition to a 7r-state around 
the region in parameter space corresponding to the Kondo regime for normal electrodes [sl, S B Q • This state is 
characterized by the presence of an unscreened spin 1/2 in the dot and by a reversal of the sign of the Josephson 
current-phase relation. This is associated with a minimum in the ground state energy at = 0^ — (j)R — tt. Outside 
the Kondo region, the ground state of the system for superconducting leads is non-degenerate with a minimum in the 
energy at (/> = 0. This state is usually referred to as the phase. The transition between the and the tt phases takes 
place through a sequence of two intermediate phases called 0' and tt' where both minima are present 

The exact phase diagram in the e — U space obtained in the ZBWM is shown in Fig. [TJ For comparison we have 
also calculated the phase diagram in the ZBWM using both the Hartee-Fock (HF) and the finite U slave boson mean 
field (SBMF) approximations [111, [l^l . The dashed lines in Fig. [T] indicate the boundary of the tt- region within these 
two approximations. One can notice that the HF approximation tends to overestimate the stability of the pi-phase for 
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FIG. 1: (Color online) e — U phase diagram of the single quantum dot case in the ZBWM. 0, 0' vr' and vr regions are indicated 
by blue, light blue, yellow and red respectively. Dashed lines indicate the boundary of the vr region when calculated using 
the HF and the SBMF approximations. The inset show these boundaries when taking into account the full spectrum of the 
superconducting electrodes with T = t. The parameters are taken in units of t = A 



lower U values while the SBMF tends to do the opposite. On the other hand both approaches yield phase diagrams 
which are in qualitative agreement with the exact result in the ZBWM. 

The fact that the exact phase boundaries in the ZBWM lie in between this two different mean field approximations 
allows one to extrapolate some conclusions to the case of superconducting electrodes with their full spectrum. Although 
the full model have been studied using more accurate methods like numerical renormalization group @ and quantum 
Monte Carlo P, a description of the complete phase diagram using these type of methods is not available. However, 
the results for the full model that can be readily obtained using HF and SBMF yield a rather precise idea of the form 
of this diagram. In fact the results exhibited in the inset in Fig. [1] for both approximations show that for the full 
model the diagram is quite similar to the the one obtained within the ZBWM when making the correspondence i — > F, 
where F is the standard tunneling rate to the leads. As in the discrete case one would expect that the boundaries 
of the exact diagram should lie somewhere in between. We therefore conclude that the ZBWM with the prescription 
i — s- F provides a good qualitative description of the ground state properties in this kind of systems. 

Finally, a word of caution should be said regarding the applicability of the ZBWM to the case of normal leads. 
While in the superconducting case the ground state properties are rather insensitive to the detailed structure of the 
leads, this is no longer the case when the superconducting gap goes to zero and the leads become normal. One can 
nevertheless rely on the ZBWM for a rough qualitative description of the system in the normal state [23j . 



In what follows we analyze the evolution of the ground state properties of systems containing an increasing number 
of sites in the Hubbard chain and within the ZBWM for describing the leads. The overall properties of these systems 
are best illustrated by the evolution of the corresponding e — U phase diagrams with increasing N . In Fig. [5] we 
show these diagrams for N — 2,3,4,5. As can be observed all these diagrams keep a certain resemblance with the 
single dot one (Fig. [1]) in the sense that there appear lines separating the zero regions for the empty (e > 0) and fully 
occupied (eo + C/ < 0) dots limit. However, one can notice a strong even-odd effect which manifests in the presence or 
absence of a tt region around the line e « — C//2 corresponding to the half-filled case. Moreover, a distinct feature is 
the appearance of new branches of vr regions with increasing N around the lines e = and t + U = Q. For even N the 
absence of a 7r-phase around the half-filled region is explained by the fact that in this case there are N electrons in 



IV. PHASE DIAGRAMS FOR THE MULTI-DOT CASE 
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FIG. 2: (Color online) e — U phase diagrams for A'^ = 2, 3, 4, 5. The color convention is the same as in Figl. 

the dots region which couple antiferromagnetically forming a singlet state with total spin S = 0. On the contrary, for 
an odd number of sites at half filling there is always an unpaired electron which leads to a ground state with S = 1/2, 
corresponding to the 7r-state. 

The other relevant feature, namely the appearance of new branches of 7r-regions around the lines e = and e + U = 
can be readily associated with the condition of a total charge in the chain region corresponding to an odd integer. 
This can be checked by analyzing the behavior of the total mean charge per spin (iVe|, Nei) of chains with = 3, 4 
as a function of the level position for a given value of U (Fig. [3|). As can be observed the 7r-regions are associated 
with the magnetic symmetry breaking {Ne-\ ^ Nei) which occurs around the regions where the total mean charge in 
the chain is an odd integer. 



V. SPIN-SPIN CORRELATIONS AND SPATIAL DISTRIBUTION OF SPIN EXCESS 

An additional insight into the nature of the correlated ground state can be obtained from the analysis of spin 
correlations. For this purpose we calculate both the intra chain spin correlations (< SiSj >, i,j = 1,N) and the 
correlation between spins in the chain and one of the leads < SLSi >. In Fig. [4] we show the average correlation 
between neighboring sites inside the chain, << SiSi+i »= J2i!^^ < SiSi+i > /{N — 1) ior N = 3,4 both in 
the normal and in the superconducting state as a function of the level position for a fixed value of U. As could be 
expected this correlation acquires an antiferromagnetic character reaching its maximum (negative) value around the 
condition of half-filling. On the other hand one can clearly notice that the effect of the superconducting electrodes 
is to increase the antiferromagnetic correlations with respect to the normal state [llj. The correlation between the 
spins in the leads and in the chain is revealed by analyzing the quantity < SlSi > which is shown in Fig. [5] for 
the same set of parameters as in the previous figure. Again in this case these correlations have an antiferromagnetic 
character reaching in the normal state their maximum negative value around the condition of half-filling. However, 
in contrast with the internal spin correlations, they are strongly suppressed around half-filling for superconducting 
electrodes due to the competing effect of pairing correlations. 

It is also instructive to analyze the spatial distribution of the spin correlations between neighboring sites along 
the whole system. This is illustrated in Fig. [5] for iV = 6 and 7 in the half-filled case. One can notice that the 
overall antiferromagnetic correlations exhibit a superimposed oscillatory behavior which can be related to the Fermi 
wavelength of the half-filled case. The suppression of the antiferromagnetic correlations in the superconducting case 
introduces important differences in the oscillatory pattern in comparison to the normal case. Due to these different 
boundary conditions the patterns have an opposite phase in one case with respect to the other. On the other hand. 
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FIG. 3: (Color online) Total charge per spin in the dots region for U — 15 for = 3 (upper panel) and N = 4 (lower panel). 
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FIG. 4: (Color online) Average first-neighbors spin-spin correlation inside the chain for the same parameters as in Fig. [3] The 
red solid (blue dashed) line corresponds to the superconducting (normal) case. 



there are some differences when comparing even and odd N. In the odd case the ideal oscillatory pattern is frustrated 
and exhibits a jump of tt in the middle of the chain. 

An interesting effect that appears in the odd N case, when the ground state is doubly degenerate with S = 1/2, 
is the non-uniform spatial distribution of the excess spin. This can be appreciated in the inset of Fig. [6] where the 
mean charge per spin along the system is shown. Again the difference between the normal and the superconducting 
case corresponds to a different boundary condition at both ends of the chain: while in the superconducting case the 
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FIG. 5: (Color online) Spin correlation between the left lead and the first site of the chain. Same parameters and color 
convention as in Fig. |4l 
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FIG. 6; (Color online) Spatial distribution of first neighbors spin correlation in the half-filled case for U — 20. The inset in the 
lower panel represents the charge per spin on each site of the system. Site number 1 indicates left electrode. We use the color 
convention of Fig. |3] 



spin excess exhibits a rather uniform oscillatory pattern along the chain, in the normal case this excess tends to 
be suppressed at both ends of the chain due to the large antiferromagnetic coupling with the leads. This state is 
analogous to what the authors of Ref. [11] refer to as two stage Kondo regime in the case of a triple quantum dot 
system. 
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FIG. 7: (Color online) Current phase characteristics for U = 10 at half-filling and A*' = 2, 4, 6 (blue, red and green lines) and 
A'^ = 3, 5, 7 (red, blue and black lines). The current is scaled by an exponential exp(Q:A'') where a has been fitted to a value 
~ 1.8. /(<;/)) is in units of eA/h. 

VI. CURRENT PHASE RELATION AND SCALING BEHAVIOR FOR LARGE N 

In this section we analyze the behavior of the current phase relation for different system sizes. The current phase 
relation is obtained from the derivative of the ground state energy with respect to the phase difference = 
{e/T%)dE{(t)) / d4). For the sake of simpUcity we will concentrate in the half-filled case. For the case of repulsive 
interactions {U > 0) analyzed so far one finds that typically the current-phase characteristic for sufficiently large 
U has a sinusoidal form ~ IcSincf) with a critical value Ic decreasing exponentially with increasing N. This 
behavior is illustrated in Fig. [71 Notice that the results shown in this figure have been scaled with an exponential 
factor exTp{aN) which nearly merges them into two universal curves one for each TV parity. We have checked that the 
value of a depends on the interaction strength U roughly linearly. These results are consistent with the prediction 
of stating that the fixed point in a renormalization group analysis for repulsive interactions and in the limit of 
an infinitely long chain, corresponds to the absence of Josephson coupling. As shown in Fig. [7] this limiting behavior 
would be reached in an oscillatory way due to the appearance of 7r-junction behavior for odd N. 




FIG. 8; (Color online) Current phase characteristics for U = —10 at half-filling and A'^ — 1,2,3,4,5,6,7 from bottom to top. 
Dcished lines correspond to even A''. The current is scaled by a factor A''. 
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In order to make contact with other predictions of field theoretical calculations for a Luttinger liquid connected 
to superconducting leads we briefly analyze the behavior of the Josephson current in the presence of an attractive 
interaction ({/ < 0). Let us remark that in this case the only stable ground state corresponds to total spin S = 0. As 
shown in Fig. [5] the current phase relation for increasing N approaches a sawtooth behavior with a critical current 
scaling as The critical current is also inversely proportional to the absolute value of the local interaction \U\. 

These results are in agreement with the predictions of Ref. [l^, [l^ [2Q| . In particular Affleck et al. have shown that 
the renormalization group fixed point corresponds in this case to the condition of perfect Andreev reflection at the 
interface. Notice that this condition is achieved in the ZBWM by setting i = A and in the half-filled case. 

VII. CONCLUDING REMARKS 

In the present work we have analyzed the ground state properties of a quantum dot array represented by a finite 
Hubbard chain connected to superconducting electrodes. We have shown that the ZBWM provides a suitable low 
energy theory for a qualitative description of these systems. This method has allowed us to obtain the phase diagrams 
in the e — U plane up to A^ < 10 sites in the array by numerical diagonalization. We have shown that the appearance 
of TT-j unction behavior requires a total odd number of electrons in the array. We have presented an analysis of the spin 
correlations and their spatial distribution which provides a more detailed insight into the ground state properties. 
Typically the presence of superconducting electrodes tends to suppress the antiferromagnetic correlations between 
the chain and the leads. While for the single dot case this suppression is associated with the disappearance of the 
Kondo effect, for quantum dot arrays the situation is more complex due to the interplay with the internal spin 
correlations of the chain. A remarkable feature found in this work is the nonuniform spatial distribution of the spin 
excess characteristic of the 7r-junction state. We have finally analyzed the evolution of the current-phase relation and 
its scaling with increasing N. For attractive interactions the results obtained are consistent with previous theoretical 
studies based on a Luttinger liquid description of the correlated region [l^, |2^. On the other hand, we have shown 
that in the half-filled case and for repulsive interactions the critical current decays exponentially with increasing A^ 
although it exhibits an oscillation in the sign due to the alternating occurrence of and tt behavior for even and odd 
N. 
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